Volume Elements of Monotone Metrics on the n x n Density Matrices as 
Densities-of-States for Thermodynamic Purposes. II 
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We derive explicit expressions for the volume elements of both the minimal and maximal mono- 
tone metrics over the (n^ — l)-dimensional convex set of n x n density matrices for the cases n = 3 
and 4. We make further progress for the specific n = 3 maximal-monotone case, by taking the limit 
of a certain ratio of integration results, obtained using an orthogonal set of eight coordinates. By 
doing so, we find remarkably simple marginal probability distributions based on the corresponding 
volume element, which we then use for thermodynamic purposes. We, thus, find a spin-1 analogue 
of the Langevin function. In the fully general n = A situation, however, we are impeded in making 
similar progress by the inability to diagonalize a 3 x 3 Ifermitian matrix and thereby obtain an 
orthogonal set of coordinates to use in the requisite integrations. 
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I. INTRODUCTION 



In this communication, we report a number of results pertaining to a certain quantum-theoretic model of the 
thermodynamic properties of a system comprised of spin-1 particles, in particular, a small number of them. Our 
approach can be contrasted with the standard (Jaynesian [^]-|3|) one, which, "in some respects can be viewed as 
semiclassical [and] can presumably be justified when the number of the constituent particles is large — in which case 
the random phases can be averaged over . . . However, in the case of a small system . . . there seems to be no a priori 



1 



reason for adopting the conventional mixed state approach" Q. Vigorous criticisms of the "orthodox information- 
theoretic foundations of quantum statistics" have been expressed by Band and Park in an extended series of paper 
1^ . Park H] himself later wrote that "the details of quantum thermodynamics are presently unknown" and "perhaps 
there is more to the concept of thermodynamic equilibrium than can be captured in the canonical density operator 
itself." Additionally, Lavenda argued (as detailed in sec. [ID) that there are deficiencies — from a probabilistic 



viewpoint — with the "Brillouin function" (used in the study of ferromagnetism) , which is yielded by the standard 
methodology Q. Contrastingly, Lavenda asserts that the "Langevin function" is free from such defects. (Aharoni ||, 
pp. 83,97,98], citing Yatsuya ei a/ in support, reaches similar conclusions under an assumption of complete spatial 
isotropy and arbitrariness of the direction of the applied field.) We have previously found App. I] (cf. [Q) that for 
the spin-^ systems, the form of analysis we will pursue here for the spin-1 systems, does, in fact, yield the Langevin 
function. Thus, our methodology appears to be not subject to the criticism of Lavenda. One of the principal results 
below will be a spin-1 version (Fig. ^ of the (spin-i) Langevin function. One possible application of these results is 
to quantum chromodynamics, where one can regard the antiscreening of the Yang-Mills vacuum as paramagnetism 
for the gluons, which are charged particles of spin 1 [[ll|Jl2| 




In the context of entangled quantum systems — to which we seek to apply our analytical approach in sec. [II — it 
has been argued that "there are situations where the Jaynes principle fails" . "The difficulties in understanding 
of the Jaynes inference scheme are due to the fact that the latter is just a principle and it was not derived within 
the quantum formalism" |p^ . Friedman and Shimony |l6|| , in a classical rather than quantum context, claim to have 
"exhibited an anomaly in Jaynes' maximum entropy prescription," cf. [p^-|l9|. 

II. SPIN-1 SYSTEMS 
A. Background 

Bloore had studied the geometrical structure of the eight-dimensional convex set of spin-1 density matrices, 
which he denoted, 

/a h g\ 

a,5,ceR, /,g,/ieC (1) 

(we incorporate the notation of ||15[ into ours). In the (MATHEMATICA [^) computations upon which we rely, a 
sequence of transformations — suggested by this work of Bloore — is implemented, leading to the full separation of 
the transformed variables In particular, we make use of a four-dimensional version {r, 01,02, S3) of spheroidal 

coordinates [p2| . 

In recent years, there have been several studies [2^-|2^ concerning the assignment of certain natural Riemannian 
metrics to sets of density matrices, such as (0). These various metrics can all be considered to be particular forms 
of monotone metrics p^-[30|]. (Contrastingly, there is only classically, as shown by Chentsov | |3l[ |, a single monotone 
metric — the one associated with the Fisher information. Morozova and Chentsov | |3^ sought to extend this work 
to the quantum domain, while Petz and Sudar [ p8| further developed the line of analysis.) Of particular interest are 
the maximal monotone metric (of the left logarithmic derivative) — for which the reciprocal (^^^) of the "Morozova- 
Chentsov function" is the arithmetic mean — and the minimal (Bures) monotone metric (of the symmetric 
logarithmic derivative) — for which the reciprocal (^^) of the Morozova-Chentsov function is the harmonic mean. 

Also of considerable interest is the Kubo-Mori metric — the reciprocal ( ^ll^^ ) of ttic Morozova-Chentsov 
function of which is the logarthmic mean - but we do not study it here, if for no other reason than that it appears to 
be computationally intractable for our purposes. In any case, we should note that in [ p5[ , relying upon a variant of 
a statistical test devised by Clarke p^ ], the Kubo-Mori metric was found to give rise to a less noninformative prior 
distribution than the maximal monotone metric, which itself was shown to be most noninformative (but only if — 
due to nonnormalizability — the pure and some collection of nearly pure states, were eliminated from consideration). 
We will make use of the Morozova-Chentsov functions in deriving (by adopting certain work of Dittmann ]3^) our 
formulas for the volume elements of the minimal and monotone metrics. 
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B. Principal Results 



The first set of questions which we wish to address is whether one can normalize the volume elements of these 
two (minimal and maximal) metrics over the ei^/ii-dimensional convex set of spin-1 density matrices, so as to obtain 
(prior) probability distributions. In the maximal case (as also in its spin-1/2 counterpart eqs. (43)-(46)]), the 
answer is strictly no, since we encounter divergent integrals. Nevertheless — similarly to the analysis in the spin-1/2 
case — we are able to normalize the volume element of the maximal monotone metric over a subset of the entire 
convex set, omitting the pure and some collection of nearly pure states. Then, by taking the limit (in which the 
subset approaches the complete set) of a certain ratio, obtain a /ower-dimensional marginal probability distribution. 
In this manner, we have been able to assign the probability distribution (cf. (fij 



15(1 -q)^A 

to the two-dimensional simplex spanned by the diagonal entries of (|l|). Of course, by the trace condition on density 
matrices, we have that a + b + c = 1. The asymmetry under exchange that is evident between a and either & or c — 
but lacking between b and c themselves — in (||), is attributable to the specific sequence of transformations, suggested 
by the work of Bloore ||l^, employed below, following his notational and analytical scheme. Consequently, it is quite 
natural to intrepret the variable a (despite the particular ordering of rows and columns used by Bloore in (|l|) , whose 
notation we have adopted from the outset of our analysis) with the middle level of the three-level system (the one 
inaccessible to a spin-1 photon, due to its masslessness), and b and c with the other two (accessible to a two-level 
system). With another, but equally valid sequence of transformations, we could have interchanged the roles of these 
variables in (^). (Let us also note that without the factor, (1 — a), (^) would be proportional to a member of the 
Dirichlet family of distributions p6[.) 

The univariate marginal distributions of (^ are (Fig. |I]) (a member of the family of beta distributions). 



15(1 -a)V^ 
4 ' 

having a peak at a = |, at which the probability density equals « 1.44338 and (Fig. |2|), 

15(l-6)(l + 35) 
32 Vb ' 



(3) 



(4) 



(and similarly for the diagonal entry c). The expected values for these distributions are (a) ~ |, and (b) = (c) = |, 

so (a) + (b) + (c) = 1. Also, (ab) - ^, (a^) = A, (^2) ^ ^^2^ ^ i_ 
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probability density 




FIG. 2. Univariate marginal probability distribution over the diagonal entry b 



C. Thermodynamic Analyses Based on Two Diagonal Hamiltonians 



1. The first diagonal Case (Xs) 



Let us now consider the observable, one of a standard set (but, due to the particular sequence of transformations 
we employ, suggested by the work of Bloore p5[ , we take the liberty of harmlessly permuting the first and third rows 



and columns of the usual form of presentation) of eight Hermitian generators of SU{3) p7 38|, 




(5) 

which might possibly function as the Hamiltonian of the spin-1 system. The expected value of (|^) with respect to (|l|) 
is (Ag) = Tr(/9Ag) = . Multiplying the univariate marginal probability distribution by the Boltzmann factor 
exp(— /^(As)) and integrating over a from to 1, we obtain the ("weak equilibrium" [||) partition function (cf. js^), 

Q{f3) = Y^(5e-^(6y^(e^'3 _ 3i(2;3+ V3)V^erfi(3^y^))), (6) 

where erfi represents the imaginary error function £I^H)_ (The error function erf(z) is the integral of the Gaussian 
distribution, that is, e~* dt.) In "strong equilibrium," zero prior probability is assigned to those density matrices 

which do not commute with the Hamiltonian, while in "weak equilibrium," this requirement is not imposed B . 

In the conventional manner of thermodynamics, we compute the expected value (((Ag))) of (Ag) as — ■ The 

result is plotted in Fig. ^. 
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FIG. 3. Expected value of (As) as a function of the inverse temperature parameter fi 

The expected value for /? = 0, corresponding to infinite temperature, is (As) = ~ —.11547. "Two physical 

conditions must be met in order for negative temperatures to arise: the subsystem possessing the negative temperature 
must be well insulated thermally from the other modes of energy storage of the complete system, and the subsystem 
energy levels must be bounded from above and below. Thus a normally populated state with probabilities proportional 
to e"^*^' can be inverted by reversal of the order of the energy levels while populations remain intact because there 
is no convenient energy sink available. Examination of the entropy gives further insight to the idea of negative 
temperatures . . . When the energy levels are bounded from above as well as below, however, zero entropy can occur for 
both minimum and maximum energy . . . No particular difficulties arise in the logical structure of statistical mechanics 
or of thermodynamics as a consequence of these negative-temperature systems" pp. 131, 132]. 

Besides (|^), the only diagonal member of the standard set of eight Hermitian generators of S'[/(3) is 




(7) 



where is used to denote a particular angular momentum matrix p6| , p. 38]. This can be viewed as the spin-1 
analogue of the Pauli matrix, 



0-3 = 2,1 



1 

-1 



(8) 



Use of this observable (tJa) in conjunction with the volume element (inversely proportional to the |-power of the 
determinant of the 2x2 density matrix) of the maximal monotone metric for spin-^ systems, has led to the partition 
function, 



Q{(i) 



sinh/3 



(9) 



yielding as the expected value, the negative of the Langevin function 

d\ogQ{p) _ 1 



d(3 



coth/?. 



(10) 



2. The second diagonal case (\z ) 

Now, we have additionally that (A3) = Tr(pA3) — b — c. Then, using a combination of symbolic and numerical 
integration, we have conducted a thermodynamic analysis for A3, analogous to that above ((||) and Fig. 0) for As. In 
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Fig. ^, we display the expected value (((A3))) of (A3), obtained numerically, along with that predicted by the negative 
of the Langevin function (10). The two curves are rather similar in nature, with the negative of the Langevin function 
assuming, in general, greater absolute expected values for a given f3. 
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FIG. 4. Expected value of (A3) as a function of the inverse temperature parameter (3, along with the expected value predicted 
by the (steeper-at-the-origin) negative of the Langevin function ( p^ ) 

D. Comparative properties of the Brillouin and Langevin functions 

1. Critique of Lavenda 

The usual use of the Langevin function is to describe the thermodynamic behavior of nonintcracting particles with 
nonquantized spin. In this same paradigm, the Brillouin function, that is tanh/3, is employed for spin-^ particles 
(Its spin-1 counterpart is i^^cosh/j ") However, Lavenda pp. 193] has argued that, in contrast to the Langevin 
function, the Brillouin function lacks a proper probabilistic foundation, since its generating function cannot be derived 
as the Laplace transform of a prior probability density". He also writes 0, p. 198]: "Even in this simple case of the 
Langevin function ... we have witnessed a transition from a statistics dictated by the central-limit theorem, at weak- 
fields, to one governed by extreme-value distributions, at strong-fields. Such richness is not possessed by the Brillouin 
function, for although it is almost identical to the Langevin function in the weak- field limit, the Brillouin function 
becomes independent of the field in the strong-field limit. In the latter limit, it would imply complete saturation 
which does not lead to any probability distribution. This is yet another inadequacy of modeling ferromagnetism by 
a Brillouin function, in the mean field approximation." Additionally, Lavenda asserts [Q, p. 20] that the Langevin 
function "has empirically been targeted as providing a good description of hysteresis curves in ferromagnetic materials 
when the field due to interdomain coupling is added to the actual internal field to produce an effective field. This 
effective field is analogous to the Weiss mean field experienced by individual magnetic moments within a single domain. 
For a sufficiently large interdomain coupling parameter, an elementary form of hystersis loop has been observed in 
the modified Langevin function. Moreover, since the generating function must be expresed as a Laplace transform 
of a prior distribution, in order to make physical as well as statistical sense, this rules out certain other candidates 
like the Brillouin function. For particles of spin-^, the generating function would be proportional to the hyperbolic 
cosine, and the hyperbolic cosine cannot be expressed as a Laplace transform of a prior distribution." 

2. The Jiles-Atherton theory 

In his claims that the Langevin function provides a good description of certain magnetic phenomena, Lavenda refers 
to the work of Jiles and Atherton . "The Jiles-Atherton theory is based on considerations of the dependence of 

energy dissipation within a magnetic material resulting from changes in its magnetization. The algorithm based on 
the theory yields five computed model parameters, Ms, a, a, k and c, which represent the saturation magnetization, 
the effective domain density, the mean exchange coupling between the effective domains, the flexibility of domain 
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walls and energy-dissipative features in the microstructure, respectively . . . The model parameter a is derived from 
an analogy to the Langevin expression for the anhysteretic magnetization Man as a function of both temperature T 
and field H for a paramagnet . . . However, in the Jiles-Atherton theory, the spin entity (m) is not an atomic magnetic 
moment m = n/Uf,, where /ib is the Bohr magneton, as in the original Langevin expression. Rather, it represents 
the moment from a mesoscopic collections of spins that we refer to as an 'effective domain;' each 'effective domain' 
possesses a collective magnetic moment (m). These effective domain entities may or may not correspond to actual 
magnetic domains" [ pO[ . 

3. Critique of Brody and Hughston 

Brody and Hughston P,|5l[ propose the use of the negative of the Langevin function, that is (|lO|), for the internal 
energy of a (small) system of spin-^ particles in thermal equilibrium. Brody (personal communication) has suggested 
that "after all, standard (Einstein's) approach to quantum statistical mechanics does seem to work for bulk objects, 
so there seems to be some kind of 'transition' from micro to macro scales" . "We note that in the case of the quantum 
canonical ensemble the heat capacity for this [spin-i] system is nonvanishing at zero temperature. Since it is known 
in the case of many bulk substances that the heat capacity vanishes as zero temperature is approached, it would be 
interesting to enquire if a single electron possesses a different behaviour, as indicated by our results" 

4. Critique of Aharoni 

In sec. 5.2, entitled "Superparamagnetism" (that is, the "phenomenon of the loss of ferromagnetism in small 
particles") of his recent text H, pp. 97, 98], Aharoni writes: "A single particle of such a small size cannot be made 
or handled. Experiments are therefore carried out on an ensemble of particles, which in most cases have a wide 
distribution of particle sizes. Such particles would give rise to a superposition of Langevin functions with different 
values of = MsV in the argument, and the measured curve could not possibly look like the Langevin function 
. . . With improved techniques for producing very small particles, their size distribution has become narrow enough for 
a pure Langevin function [citing |^] to be observed . . . The calculation [of the Langevin function] can now be said to 
have been confirmed by direct experiment. Of course, a Langevin function (or any other similar function) can always 
be fitted to such data for a rather narrow temperature range [citing [ p| , p2[ ], but the remarkably narrow distribution 
of [ [^] can be fitted to such a function over a wide temperature range. In this respect this experiment is still quite 
unique in the literature" (cf. [^^). 

Earlier ||, p. 84], Aharoni asserts: "And there is no mistake in this algebra: there are only two differences between 
this calculation [of the Langevin function] and the study of a gas of paramagnetic atoms in section 2.1 [yielding the 
Brillouin function]. One is that the function 6 [the angle to the fixed magnetic field] is continuous here, while this 
variable had discrete values in section 2.1 and the other is that the magnetic moment /i was that of a single atom 
there, while here it is the moment of a large number of atoms, coupled together. However, the second difference is 
only quantitative and not qualitative, and the first one should not make any difference, especially since the energy 
levels of a large spin number S are very close together, and look like a continuous variable. It is thus true that if there 
was no other energy term besides the isotropic Heisenberg Hamiltonian, it would have been impossible to measure 
any magnetism in zero applied field, and there would be no meaning to a Curie termperature, or critical exponents, 
or any of the other nice features mentioned in the previous chapters. Theorists who calculate these properties never 
pay attention to the fact that the possibility of measuring that which they calculate is only due to an extra energy 
term, which they always leave out. Of course, a magnetization as in [the Langevin function], which is zero in zero 
applied field, contradicts not only experiments ... It is also in conflict with everyday experience ... It is because real 
magnetic materials are not isotropic, and not all values of the angle 9 are equally probable." 

5. Relations to modified Bessel functions 

It is also interesting to observe that the Brillouin and Langevin functions are both instances [D — \ and D — 
respectively) of the two-spin correlation function of the _D-vector model (of isotropically-interacting D-dimensional 
classical spins) for a one-dimensional lattice [pq. Fig. 2.3], 
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Io{J) 



where I{x) is a modified Bessel function of the first kind. 



E. Derivation of the Volume Element of the Maximal Monotone Metric for Spin-1 Systems 

We now discuss the manner in which our first reported result (|^) was derived. To begin with, we noted that Bloore 
| p5[ had suggested the transformations ("to suppress the dependence on [the diagonal entries] a,b and c by scaling 
the [off-diagonal] variables f,g,K^), 

f = Vb^F, g = y^G, h = ^/abH. (12) 

The positivity conditions on the density matrix p, then, took the form ||l5| , eq. (15)], 

\F\<1, |G| < 1, \H\<1, + |G|2 + - 2Rc(FGiJ) < 1. (13) 

Bloore indicated that if = Fr + \Fi were simply real (that is, Fj = 0) , then, the last condition could be rewritten 

as 

l-Fl> Gl + G] + Hi + H] - 2Fr{GrHr - GiHj). (14) 

Then, by the rotation of the (G^, i/^J-pair of axes through an angle of j and the (G/, -ff/)-pair of axes, similarly, the 
allowed values of Gr, Hr, Gj, Hj could be seen to lie inside a four- dimensional spheroid, two of the principal axes of 
which had a length equal to \/l + Fb. and the other two, Vl — Fr (so the spheroid is neither predominantly "oblate" 
nor "prolate" in character). Bloore's argument easily extends — using two additional identical rotations, not now 
necessarily however, equal to j — to the general case, in which F is complex. (If we were to single out G, say, rather 
than F, we would obtain an analogue of (H), containing b, not a, as its distinguished variable among a, b and c.) 

We proceed onward by converting to polar coordinates, Fj — s cos v, Fr = s sin v. Then, we rotate two of the six 
pairs formed by the four axes obtained by the two ^-rotations — each new pair being comprised of one member from 
each of the two pairs resulting from the ^-rotations suggested by Bloore. The two new angles of rotation now both 
equal ^ cot~^(tani'). We have, after performing the four indicated rotations, transformed the resultant set of axes 
( Ji, J2, J3, J4), using a set of four "hyperspheroidal" coordinates (r, 01,62,63), 

Ji = rVT+scos^i, J2 = r Vl + s cos 6*2 sin di , (15) 



J3 = rVl — s cos 63 sin 62 sin Oi, J4 = r\/\ — s sin 63 sin 62 sin 6\ . 

The Jacobian of the total transformation (the scaling transformations ([T^), the four rotations of pairs of axes, along 
with the introduction of polar and hyperspheroidal coordinates) is, 

J{a, b, c, s, r, 61, 62) = a%'^c^r^s{l - s^) aindi^ sin 6*2. (16) 

In the new variables, the determinant of the 3x3 density matrix p takes the simple form (being free of the four 
angular variables — u, 61,62, 63), 

IpI =a6c(l-r2)(l-s2). (17) 

To arrive at the volume element of the maximal monotone metric, which we seek to integrate over the eight- 
dimensional convex set of spin-1 density matrices, we adopted an observation of Dittmann ]35| ] regarding the eigen- 
values of the sum (Lp + Rp) of the operators of left (ip) and right multiplication (i?p) for the minimal monotone 
metric. He had noted that for the nx n density matrices, in general, these eigenvalues would be of the form pi +pj 
{i,j = 1, . . . , n), where the pi {i = 1, . . . ,n) are the eigenvalues of p itself. Then, we observed [Q (cf. |5^, eq. (24)]) 
— making use of the fact that the determinant of a matrix is equal to the product of eigenvalues of the matrix — 
that the corresponding volume element would be proportional to the square root of the determinant of {Lp + Rp)^^ 
(or, equivalently, to the reciprocal of the square root of the determinant of Lp + Rp). We adopted this line of ar- 
gument to the maximal monotone case by replacing the sums pi + pj by (twice) the corresponding harmonic means 
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(2piPj/{pi +Pj)) 1^^. Then, the volume element of the maximal monotone metric can be shown to be proportional to 
the product of \p\~^ and a term W2 — \p\, where W2 denotes the sum of the three principal minors of order two of 
p. (The spin-^ counterpart of this product — yielded by the analogous argument, as well by an independent one |2^ , 
eq. (3.17)] — is simply For the minimal monotone metric in the spin-1 case, on the other hand, the volume 

element — now based on the terms Pi +pj — is proportional to \p\~^ {w2 — The occurrence of the term W2 — \p\ 

in the denominator of the integrand — and not the numerator, as in the spin-1 maximal monotone case — results 
in a substantially more difficult multiple integration (cf. ([l8|)), for which we have been able to exactly integrate over 
the variables u, 03,82, leaving us with an expression containing an elliptic integral of the first kind. However, further 
attempts to use numerical methods for this five-variable expression have so far not been successful. 

The spin-i form of the volume element for the minimal monotone metric is proportional to |/o|^2. Its use for 
thermodynamical purposes, led j57j to a somewhat different ratio (cf. -/2 (/3)//i (/3) , of modified Bessel functions, 
than the negative of the Langevin function (p^, which is expressible as -l3/2iP)/ Ii/2iP)- (Brody and Hughston 
I psf had, first, arrived at the former ratio, but then [^,^ [personal communication], by using "phase space volume" 
rather than "state density' [i. e., weighted volume]", concluded that the latter ratio was more appropriate. However, 
for higher dimensional cases — spin-1, the case under investigation here, being the simplest such — Brody and 
Hughston stated that "the energy surfaces are not fully ergodic with respect to the Schrodinger equation, and thus 
we cannot expect to be able to deduce the microcanonical postulate directly from the basic principles of quantum 
mechanics.") So, it can be seen that as we pass from the spin-^ case to the spin-1 situation, for both the minimal 
and maximal monotone metrics, we have to deal with an additional term (besides a certain power of |p|) of the form, 
W2 — p, to some power. 

To find the bivariate marginal probability distribution (^, we first performed the eightfold multiple integration 
(using the trace condition to set c = 1 — a — 6), 



pL pL — a po p ti pTT pTT pZTT p ZlZ 

/ / / // /// W'^{w2-\p\)J{a,h,c,s,r,ei,e2)AvAe:ide2dei<lrdsdbda. (18) 
Jo Jo JO Jo Jo Jo Jo Jo 

Then, we (following established Bayesian principles for dealing with nonnormalizablc prior distributions |59|-|6l|]) took 
the ratio of the outcome of the (initial sixfold) integration — that is, after all but the last two stages (over a and 6) 
— to the result of the complete integration. In the double limit, i? — > 1, S' — *■ 1, the ratio converges to the probability 
distribution (^) over the two-dimensional simplex. (The result was invariant when the order in which the limits were 
taken was reversed.) This computational strategy was made necessary due to the fact that the integration diverged 
if we directly used the naive upper limits, r = 1 and s = 1, evident from (|l7|). A similar approach had been required 
(due to divergence also) in the spin-i case, based on the maximal monotone metric |l^, eqs. (43)-(46)]. There, the 
corresponding volume element (cx |p|~^) was integrated over a three-dimensional ball of radius R. As i? ^ 1, this 
ball coincides with the "Bloch sphere" (unit ball) of spin-i systems. A bivariate marginal probability distribution 
was then obtainable by taking the limit i? — > 1 of a certain ratio — with the resultant univariate distribution simply 
being uniform in nature, and leading, when adopted as the density-of-states, to the Langevin model (p^. 

Since it has been argued elsewhere [|5^,Q that of all the possible monotone metrics, the maximal one is the 
most noninformative in character, it would seem plausible that the maximal metric might be singled out to provide 
density-of-states (structure) functions for thermodynamic analyses. (For additional results pertinent — in the context 
of "universal quantum coding" — to comparative properties of the maximal and monotone metrics, see ||6^.) 

We would also like to report a success in performing the eightfold multiple integration in (|l^), after reordering the 
individual integrations so that the ones over r and s — rather than u and 63 — were performed first. Then, we were 
able to compute the limit of the ratio of the result after the first two integrations to the complete eightfold integration. 
This gives us a sii-dimensional marginal probability distribution (free of r and s, and independent of the particular 
values of and ^3), 

15(1 -a)V^ sin 6>iSin 6^2^ 

One of the two-dimensional marginal probability distributions of (|l9|) — the one over the simplex spanned by the 
three diagonal entries (a, b, c) of the density matrix (|l|) — is our earlier result (|2|). 

The focus of the study here has been on the application of symbolic integration to the modeling of the thermody- 
namics properties of spin-1 systems. With the application of numerical integration methods, however, we might hope 
to broaden our work to encompass non-diagonal Hamiltonians (such as Ai, A2, A4, A5, Ag, A7). 
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F. Thermodynamic Analysis Based on a Non-Diagonal Hamiltonian (Ai) 



Let us proceed with a "strong equiHbrium" analysis Q of the non-diagonal observable, 



Ai 




(20) 



(It would appear that all our conclusions of a thermodynamic nature should be precisely the same for the other five 
non-diagonal observables, as well.) Rather than considering the full eight-dimensional convex set of 3 x 3 density 
matrices, however, we restrict our considerations to the two-dimensional convex set of 3 x 3 density matrices which 
commute with Ai — so, we are considering the case of "strong equilibrium" |p|. We, then, assign the volume element 
of the Bures (minimal monotone) metric to this set (using formula (10) in p^), as a density-of-states. (Since the set 
is composed of commuting matrices, we are in a classical situation in which there is a unique monotone metric [ p8[ , 
so there is, in fact, no distinction between maximal and minimal ones, as earlier.) Our partition function then takes 
the form, 



Q(/3) = 



"'0 



1-5 



-I3(l-2g-h) 



(We perform the inner integration symbolically, and the outer one, numerically.) In Fig. 
associated expected value, along with that given by the negative of the Langevin function 



(cf. Fie 



(21) 

we plot the 



expected values 
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FIG. 5. Expected value of (Ai) and the negative of the Langevin function ( [lO|) 

Since these two curves quite remarkably almost coincide, in Fig. || we plot the result obtained by subtracting from 
the negative of the Langevin function (|l0|), the expected value of Ai. 
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difference of expected values 
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FIG. 6. Difference between the negative of the Langevin function (|lO|) and (Ai) 



III. SPIN- 1 OR COUPLED SPIN-i SYSTEMS 

A. Required Transformations of the Fifteen Parameters 

The possibility of extending our line of analysis to the {n^ — l)-dimensional convex sets of n x n density matrices, 
for n > 3, obviously, remains to be fully investigated. The 3x3 density matrix (|l]) can be considered to be embedded 
in the upper left corner of the n x n density matrix. Then, for this block, we can employ precisely the same set of 
transformations as in the analysis above. (Of course, there are "^"g""'"-' ways in which to actually position the 3x3 
block. Depending upon which we choose, we would expect to break certain symmetries between the diagonal entries, 
as we witnessed with the probability distribution (|^).) This leaves us with n — 3 new diagonal entries and — 10 
new off-diagonal variables. 

For the case n = 4 (which we had previously studied for various scenarios in which most of the fifteen parameters 
were ab initio set to zero), 



a,b,c,deR, f,g,h,o,p,qeC (22) 



we have found that in the determinant, the cross-product between the real and complex parts of each new (scaled) 
off-diagonal entry (that is, o,p,q) is zero. Then, the extension of our earlier reasoning led us to form two triads^ one 
element of each triad coming from one of these three pairs, and attempt to rotate them (in three-dimensional space), 
so as to eliminate the non-zero cross-product terms [|3 - 66 . We have also investigated the possibility that by rotating 



/a 


h 
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c 


q 


\o 


P 


9 


d 



the six variables in question in six-dimensional space, one could, in a single process, nullify all the non-zero cross- 
product terms. To accomplish this, it would be necessary to diagonalize a 6 x 6 symmetric matrix (or, equivalently, 
a 3 X 3 Hermitian matrix). The three pairs of equal eigenvalues of this matrix are given by the solution of the cubic 
equation, 

+ ^((6 " - r^s - 2s^ - r^s{l - 4sin6'i^ sme2^))x^ + 3(1 - r^){l - s^)x + (1 - r^f{l - s'^f = 0. (23) 

There does not appear to be any very simple general form for the three roots of this equation. (In principle, we are 
able to diagonalize the 6x6 matrix, but the computations required to fully implement the associated transformations 
of the six variables, seem to be quite daunting.) However, if we set the factor (1 — 4sin0i^ sin 02^) to zero (which can 
be accomplished by taking either 6i or 62 to be zero) , then the three roots are of the form 

x= (-l + r2)(l + .s), x = i(l-s)(-2-.s+p), x = -i(l-s)(2-t-s+p), (24) 
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where 



p = _^4r2(l + s). (25) 
If we set the factor to be -3 (by taking both Oi and O2 to be -1), then the roots are 

x= (-l + r2)(l-5), = i(l + s)(-2 + .s+p), a, = -1(1 + s)(2 - s + p), (26) 

where now 



p = _^ 4^2(1 _ (27) 

The situation becomes simpler still if we set the parameter s to zero (so that the 2x2 submatrix of (|^) in the upper 
left corner is diagonal in nature). Then, the three roots are —1 — r,— 1 + r and — 1 + r^. Similarly, for r = 0, the roots 
are — 1 — s,— 1 + s and — 1 + . If we set s = 1 (corresponding to a pure state) , two of the roots are zero and the third is 
—2(1 — + sin 611^ sin 6*2^). For r = 1, there are also two zero roots and one equalling —2 + s + — 2s sin6'i'^ sin6'2'^- 
Let us proceed to analyze in detail the case r = 0, for which the three eigenvalues, as just noted, are —1 — s,— 1 + s 
and — 1 + s^. Setting r to zero is equivalent (cf. (p^)) to nullifying the entries g and h in (p^, so the scenario is 
eleven-dimensional in nature. We have found the eigenvectors corresponding to these three eigenvalues, that is, for 

cos V sin V 1 , , sin v cos v \ ^ , , 



for (-1 + s), 



and for (— 1 + s^) 



cos V sin V 1 , , sin v cos v 1 



(0,1,0,0,0,0), (1,0,0,0,0,0). 



We transformed the three complex variables o,p and g, on the basis of these eigenvectors, to a new set of six real 
variables {Ki . . . , Kq)^ between which there are no nonzero cross-product terms in the determinant of the transformed 
matrix. Then, we employed a six-dimensional version of spheroidal coordinates (w, Cii ^3, C4j Cs) (cf- (|l5|)), 

Ki — vVl + scos^i, K2 — vVl + s cos ^2 sin f 1 , (29) 



— v\/l — scos^3 sin ^2 sin^i, K4 = vy/l — s cos ^4 sin ^3 sin ^2 sin^i, 

— V cos ^5 sin ^4 sin ^3 sin ^2 sin^i, Kg = v sin ^5 sin ^4 sin ^3 sin ^2 sin^i. 

The determinant becomes, then, simply abcd{l — s'^){l — v'^). The Jacobian corresponding to the sequence of trans- 
formations is (cf. (p^)) 

J(a, 6, c,fi,s,u, ^1,^2,^37 '?4) = ab^c^d'^s{l — s^)u^ sin^i^ sin^2'^ sin^3^ sin^4. (30) 



B. Volume Elements of the Maximal and Minimal Monotone Metrics for Spin-| Systems 

The integrand based on the maximal monotone metric is (employing once again our ansatz based on the work of 
Dittmann p5[ , in which we replace the arithmetic means of the eigenvalues by their harmonic means, in accordance 
with the corresponding Morozova-Chentsov functions), the product of this Jacobian and the ratio of the term {W2W3 — 
w'^ — \p\) to , where W2 is the sum of the six principal minors of order two and is the sum of the four principal 
minors of order three. We performed the integration over the eleven-dimensional convex set, the innermost integrations 
being over v from to V, and s from to S. (The integration would diverge if we were simply to set either S or V 
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to 1.) Then, we were able to integrate over the six angular variables (u,^'s), but not fully over the three-dimensional 
simplex spanned by the diagonal entries (a, 6, c, d). 

The problem, in this regard, seems attributable to the fact that the powers of a, b and c — but not d — in the 
Jacobian (30) are all less than three. If all four powers were, in fact, equal to three, then dividing by with 
\p\ cx abed, would yield a term of the form {abcd)~2 — similar to the term {abc)~^ , encountered in the successful 
spin-1 analysis above — which would be integrable over the simplex. (We would, in fact, obtain such a desirable 
term if we could otherwise proceed with a full fifteen-dimensional analysis.) However, we have the resultant term 
a~2 2 c~ 2 (i~ 2 , which appears to be not integrable. Thus, we have not so far been able to extend the form of 
analysis taken here for the spin-1 systems to the convex set of spin-| or, equivalently, coupled spin-^ systems (cf. 
I pof ) . (Following the lead of Bloore |l^ , one might also study the presumably simpler case of the nine-dimensional 
convex set of real 4x4 density matrices.) 



C. Strong Equilibrium Thermodynamic Analysis of a Non-Diagonal Hamiltonian 

Let us, in the manner of sec. |ll F| , restrict our attention to those 4x4 density matrices which commute with a 
certain non-diagonal Hamiltonian (so we remain within the framework of "strong equilibrium" ) , which we choose here 
to be 

/O 1 

10 10 

10 1 
Vo 1 




Then, after assigning the volume element of the Bures metric to the three-dimensional convex set of the mutually com- 
muting density matrices, we have found the expected value of the observable (|3l]), as a function of the thermodynamic 
parameter /3, to be as indicated in Fig. ^. 

expected value 




b 



FIG. 7. Expected value of the non-diagonal observable (|3l| 



IV. CONCLUDING REMARKS 



In this communication, we have sought to pursue a particular chain of reasoning (which, for spin-^ systems, is 
known to yield the Langevin function, thus according fully with detailed arguments of Lavenda [Q) for spin-1 and 
spin-| systems. The fundamental problems encountered en route include the derivation of the volume elements of the 
(noninformative) maximal monotone metric(of the left logarithmic derivative) for these two scenarios. By applying 
principles used in Bayesian reasoning |59|-|6l[|, we have succeeded through extensive analysis and computations in 
obtaining a six-dimensional marginal probability distribution ( p^ from the (nonnormalizable) volume element for the 
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eight-dimensional convex set of spin-1 systems (sec. ||). We have not been quite as successful, however, in addressing 
the spin- 1 scenario, in which the corresponding convex set is fifteen-dimensional in nature (sec. [II). 

We are not able to actually demonstrate, but can only conjecture, that our results have physical implications (cf. 
Q). The most questionable and debatable assumption underlying our analysis, perhaps, is that the volume element 
of the maximal monotone metric can be employed as a density-of-states for thermodynamic purposes. 

The initial impetus to pursue this line of research was provided by an extended series of forcefully-argued papers 
of Band and Park Q. In them, for various conceptual reasons, they expressed dissatisfaction with the conventional 
( " Jaynesian" ) approach to quantum statistical thermodynamics and sought to develop a conceptually superior alter- 
native approach. We believe that our analysis is quite consistent with the spirit of their work, and serves to implement 
it for certain low-dimensional scenarios. We have also been encouraged in our work by the detailed remarks of Lavenda 
0, regarding the propriety of using the Langevin function for spin-^ systems, since our methodology, in fact, yields 
this function. Also, the analyses of Brody and Hughston p|j5l|] , though differing from ours in several technical respects 
— in particular, focusing on pure states — possess similar motivations. 

It also appears natural that, amongst the continuum of monotone metrics |2^j2^ , the maximal one should play the 
distinguished information-theoretic role we have explored for it, since its noninformative nature has been independently 
established js^js^. The fact that, as a result of lengthy and demanding computations based on the volume element of 
the maximal monotone metric, we were able to arrive at the quite simple results (H) and ( [l9| ) also assists, we believe, 
in validating the rather unconventional course pursued here. 
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